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Conventional one-dimensional oscillation theorem is found to be violated for multi-component 
Schrodinger equations in a general case while for two-component eigenstates coupled by the sign- 
constant potential operator the following statements are valid: (1) the ground state (v = 0) is not 
degenerate; and (2) the arithmetic mean of nodes m, 712 for the two-component wavefunction never 
exceeds the ordering number v of eigenstate: (rii + n2)/2 < v. 



I. INTRODUCTION 

For a large number of physical applications there are 
useful the so called oscillation theorems 0, 0] that de- 
scribe some remarkable properties of solutions of the one- 
dimensional (ID) Schrodinger equation. As for ID prob- 
lems the bonded states are nondegenerate one may enu- 
merate the wave functions of stationary states by their 
numbers v in increasing energy order (v = 0, 1, 2, • • • ). It 
is known (see Q, Ch.XIII.3 in or Ch.VI, §6 in 0) that 
the eigenfunction ip v corresponding to the v-th eigenvalue 
E v vanishes exactly v times (another words, the ip v pos- 
sesses exactly v internal nodes). These statements may 
be clarified within the quasi-classical (WKB) approxi- 
mation (see 0, [Bj]). For the many-dimensional problems 
these statements are not so simple Q. 

In present study we consider the nodal properties 
of multi-component wave functions for the Schrodinger 
equation being interesting within the physical problems, 
where adiabatic approximation shell be used in modified 
forms (see e.g. Q). For example, the spin-orbit coupling 
effects also arc resulted in a multi-component diatomic 
vibrational Schrodinger equation 0, 0] • 

Let us consider the two-component Schrodinger equa- 
tion for the 2x2 matrix Hamiltonian H: 



H(p v = 



hi Via 
V21 hi 



= E v 



(Al) 



The components hi ,2 of H in Eq. (|Alj) have the tra- 
ditional form hk = ~-^-§^2 + Uk(x) with some effec- 
tive mass \x and real- valued potentials Uk (here the ro- 
tational part of energy is included in U^). The inter- 
action operators are supposed to be some real-valued 
potentials V~2i(x) — Vi2(x), and the energy levels E v 
are enumerated in increasing order, taking degeneracy 
into account, by traditional vibrational quantum num- 
ber v = 0, 1, 2, ••• . 

For molecular problems one may suppose the compo- 
nents fi2i x ) °f a vector-function tp v to be real- valued 
functions defined for some region f2, being the same for 



both components. It is supposed here f2 = [a, /?] C M 1 
for some finite (or infinite) a and f3. When one of the 
ends of the interval [a, /J], for example, a, is finite, one 
supposes that both components of solution tp v obey the 
Dirichlet boundary condition at this point: tp\ 2 { a ) = 0. 
In any case only bound states are considered and solution 
ip v is supposed to be normalized by condition 



<v^}H^iK> + <^>2> = i- 



(A2) 



The analysis of the problem in the WKB- or/and "adi- 
abatic" approximations (of course, being modified for the 
systems considered) and by the perturbation theory (for 
small enough V12) gives a wide range of possibilities. In 
particular, it is easy to find some special situations when 
the problems under consideration differ essentially from 
the traditional one-component ones. So it seems to be 
interesting to describe some general statements, that are 
valid for a large part of practical applications. 



II. SOME EXAMPLES AND GENERAL NOTES 



In general the problems of the form (|A1IA2[) seems 
to be not essentially differs from the one-component 
Schrodinger ID-equation. It may be rewritten as a sys- 
tem of two equations with the same normalization con- 
dition (lA2l): 



hupl + V12V2 = E v <pt 



(A3) 



It is clear, that when interaction potential vanish 
(Vi2{x) = for any x) the problem (|A3j) is reduced to 
some one-component eigenvalue problems for Hamilto- 
nians ft. 1,2- In this situation the Hamiltonians hi and 
/12 with essentially different densities of the bond states 
shows it clearly, that the number of nodal points for one- 
component functions does not define the number (that 
is, the v value) of the corresponding energy level. 
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If the spectra of hi and are identical and V\i = 0, 
any energy level is degenerate one, including the ground 
state level. One may also give a simple example of the 
situation, when degeneracy of all the states in the two- 
component problem Eq. flAip appears for nontrivial Vi2- 
Really, when hi = h,2 = h is the Harmonic Oscillator 
Hamiltonian on the R 1 axis and Vi2(x) = x, it easy to 
transform Eq. (|Aip to a system of two eigenvalue prob- 
lems for lD-Hamiltonians h + V12 and h — V12, that de- 
scribe the two equivalent oscillators with shifted centers. 
Hence, all the energy levels of the problem considered are 
doubly degenerate ones. Note also, that the small per- 
turbation of the problem above allows one to construct 
the two-component problems, where the components of 
the ground state have nodes. 

The simple examples above described demonstrate it 
clearly, that the nodal structure of the many-component 
wavefunctions is a much more difficult problem, than 
that one for one-component ID-system. Nevertheless, 
one may note, that in all the examples presented the in- 
teraction potentials V12 vanish, at least, in some points. 
If one exclude the situations, when V12 change its sign 
or vanishes, that is, when the interaction potential V12 is 
invertible, the situation changes drastically and one may 
prove some interesting statements on nodal structure of 
the multi-component wavefunctions. 

For invertible potential Vi2(x) one may use the sys- 
tem Eq. (IA3j) to find ip\ as a solution of the 4-th order 
differential equation (similar relation holds for tp^): 



[(h 2 - E°)V{?{hi -E v ) + V 12 M = 0. (A4) 

In particular, the most of statements below use the fol- 
lowing idea, being the result of the well known property 
of the Cauchy problem: 

Lemma. For invertible potential Vi2{x) the compo- 
nents of any nontrivial solution ip v of the problem (jAlj) 
cannot be identically zero into some open subregion of 

n. 

For the ID one-component problems it is impossible 
that the wave function and its gradient vanish simul- 
taneously at some point. In particular, near the nodal 
points the wavefunction changes the sign. For the multi- 
component problems situation is not so simple. In prin- 
ciple, for solutions of Eq. (IA4[) it is possible that at some 
point some component of the wave function only touches 
the axis. When the interaction potential does not vanish 
at this point, it follows from Eq. (|A3| that both com- 
ponents tfi 2 (x) vanish along with their second deriva- 
tives. This special situation is not stable with respect 
to small perturbations and we ignore it here. That is, 
we consider as the nodal points only the ones for which 
functions change their signs. Similar issues for one com- 
ponent problems for £1 = R"(n > 1) enforce to analyze 
the regions of a constant sign for wave functions [2[ . 
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III. USEFUL STATEMENTS 

Here we prove three useful statements that gives the 
background for discussions of the nodal properties of 
wavefunctions in the next Section. The notations in- 
troduced in this Section are used here and in the next 
Section without special comments. In particular, we use 
the symbol [n/2] for the integer part of the n/2 value. 
The following relation is important for our discussions in 
this text: 



(A5) 



An other construction is typical for the problem under 
consideration (see, e.g. [2j). Let <p v be an eigenvector of 
the two-component matrix Hamiltonian H with energy 
E v and components ip\ and ip% • Let us analyze the re- 
gions of a constant sign for the function ip\ (or for tp^), 
that is, the intervals on R 1 where the function has a def- 
inite sign and changes it at the ends. The number Kj of 
such regions corresponds to rij — Kj — 1 nodal points of 
^(j = l,2). 

Let us suppose that ip\ has Ki regions Q,j (j — 
1,2,..., Ki) of a constant sign. One may define for each 
of regions fi^ the function ip^ that equals to <p\ on 

and equals to zero out of , Let us suppose that Lp v 2 
(2) 

has K2 regions il k (k = 1, 2, . . . , K2) of a constant sign 

and (p k are the functions similar to the ones for the first 
component of ip v . Hence the two-component function ip v 
gives rise to construct the linear space K of {Ki + K 2 ) 
-dimensional vectors C defined by any set of Ki coeffi- 
cients Cj and set of K2 coefficients du and corresponding 
vector-functions x with components of the form 



-R"i 



3=1 



(i) 



X2 



K 2 

E 

fc=i 



^ (2)- 



(A6) 



In particular, \ — f v when d^ 



1 for any k and 



j values. Note also, that X1.2 and their gradients have 
to be the square integrable functions on D, as Lp\ 2 fulfill 
these requirements (the sets of a zero measure have no 
importance here). 

Now we consider the properties of the matrix M, de- 
fined for the above mentioned space K, by the vector- 
function x i n the following way: 



(C,MC) 



E 

3\k 



(c 3 -d k ) 2 {^\Vi2\<p k 2) ). (A7) 



Our attention to the matrix M is explained by the 
following reason. The further discussion is based on the 
simple expression of energy functional E{x) associated 
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with the Hamiltonian H, for the vector-function \ with 
components of the form (|A6I) . 

Statement 1. If ip v is a solution of Eg. (IAl[) . then for 
any vector CfiK one may write 



E(x)=E v (x\x)-(C,MC) 



(A8) 



Proof. The use of a standard form for the energy func- 
tional (with kinetic energy expressed as squared gradi- 
ent of wavefunction) and integration by parts enables to 
write, due to Eq. (|A6|) . the following relation: 



E(X) 



j=i k=i 



+2 



(2), 



where integration in each of integral in the sums is done 
over the regions fi^ , f2^ or Q^|~)f2j[. , respectively. 
The use of Eq. (|A3|) gives immediately 



Proof. According to definition (see Eq. (IA7[) ). the 
quadratic form (C, MC) is determined by a sum of matri- 
ces of rank 1, each of which has only one nonzero eigen- 
value 2{<pf ) \V 12 \^ ) ), if any. According to the State- 
ment 2, among those eigenvalues there are no more than 
(K1+K2 — 1) non-zero values. A half of them is negative, 
while another half is positive. In any case, according to 
the Eq. (|A5[) the largest number of strictly negative inte- 



grals is no more than [ ^j^ 2 ] . 

According to the minimax principle (see, e.g. Ch.VI. 
of [1] or Sect. XIII. 1 of Q), this means that no more than 
[ Kl t K2 ] eigenvalues of matrix M are strictly negative. 
Hence, the remaining eigenvalues of M are non-negative 
and their number is no less than 



~K X + K 2 ~ 






2 




2 



(Ki + K 2 ) - 



as follows from Eq. (|A5 



IV. THE NODAL POINTS IN 
TWO-COMPONENT PROBLEM 



(All) 



k=l 



3,k 



+2EcA<^-Virf ) ) 

(A10) 



The use of definition Eq. (|A7|) shows it clearly, that this 
is exactly expression (|A8j) . ■ 

It is interesting, that the number of positive/negative 
eigenvalues of the matrix M is closely connected with 
numbers 711,2 of nodes for components of the vector- 
function (p v . To prove it one shell note the following 
simple statement. 

Statement 2. There are no more than (Ki + K 2 — 1) 
regions f] fl^ C fl, where the product ip\tp 2 con- 
serves its sign. 

Proof. In general there are no more than n\ + n 2 inter- 
nal points of the region fl where <f\<f 2 changes the sign. 
As rij = Kj — 1, there are no more than {K\ — 1)+(K 2 — 1) 
nodes for ifi^p 2 , and, hence, there is not more then 
{K x - 1) + (K 2 - 1) + 1 = (K 1 + K 2 - 1) non-empty 
regions of the type fij fl^fc where the function ip1ip 2 
conserves the sign. ■ 

The most important result of this Section is the fol- 
lowing statement. 

Statement 3. If the sign of the potential ^12(3;) is 
fixed, then the matrix M defined by Eq. (|A7p has no less 
than \ Ei±^2±l j non-negative eigenvalues. 



For ID problems the wavefunction of the state with 
vibrational number v = 0, 1, 2, . . . has exactly v nodes 0, 
SH. Whendim(fi) > 1 , for the one-component problem 
one may only prove that, for a given state, the number of 
regions of a definite sign is not too large [2j . For the two- 
component system considered here one may prove some 
analogue of these statements. The proof is similar to 
the one in the one-component case, but it requires some 
special modifications. 

Theorem. For the v-th stationary state of the prob- 
lem Eq. (|Aip with the interaction potential V\ 2 (x) of a 
definite sign the number of points n\ and n 2 , where the 
components tp\ and Lp^ of the vector eigenfunction ip v 
change their signs, the following inequality holds: 



n\ +n 2 + l 



(A12) 



and, hence, 



ni + n 2 



< v. 



(A13) 



Proof. Let us consider the space /C associated with 
regions of constant sign for components of the wavefunc- 
tion ip v with energy E v , as was described in Section [TT1 
Each vector C € K. may define some vector-function x 
(see Eq. (IA6|) ). Due to Statement 1 the corresponding 
energy functional value can be written as 



^(X) = ^{xlx)-(C,MC). 



(A14) 
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According to the Statement 3 (remember, that Kj 
rij + 1), the matrix M has no less than 



T : 



K x + K 2 - 


-1" 




ni + n 2 + 3 


2 




2 



n\ + n 2 + 1 



1 (A15) 



non-negative eigenvalues, each of them being asso- 
ciated with some vector C and corresponding vector- 
function x- 

Let us suppose T > v + 2. Then one may find at /east 
two independent normalized vectors X^ 1 (or, equiva- 
lently, C^ 1 ' 2 * 1 6 /C being orthogonal to v exact solutions 
<fP, cp , ... , <^" _1 of Eq JAlj) with lowest energies. 

For these vectors (C®,MC®) > (J = 1,2) and, 
hence, -E(x^ 1:2 ') < E v . It follows from the variational 
principle (see, e.g. Ch.VI. of 0] or Sect.XIII.l of Q) 
that both functions jc an d x^ are solutions of Eq. ljAlj) 
with energy E v and one may construct, as a linear com- 
bination of these vector-functions, a non-zero solution x 
of the Schrodinger equation, for which, at least for one 
of regions Qj or the corresponding component of 
the solution (xi or X2) identically equals to zero. This is 
however impossible, according to the Lemma of Section 
|T1 Hence, T < v + 2, or, equivalently, T < v + 1. Hence, 
one may write 



m + n 2 



< 



< v. 



(A16) 



due to Eq. (|A15[ ). The usage of Eq. (IA5[) proves the state- 
ment. ■ 

It is clear, that (ni + n 2 )/2 is the arithmetic mean of 
number of nodes for components of the two-component 
solution. Hence, the Theorem statement Eq. (|A13[ ) is an 
analogue to the usual one-component statement for the 
case of the Schrodinger Equation in W n (n > 1) (see Ref. 

0)- 

Corollary 1. For the Theorem's conditions, the com- 
ponents of the function tp" =0 conserve their signs into 

n. 

[Indeed, for v — inequality Eq. (IA12p means n\ — 
n 2 = 0]. 

Corollary 2. For the Theorem's conditions for the 
ground state and 1^12(2;) < the signs of tp\ and ip 2 
coincide, while for V± 2 (x) > the signs mentioned arc 
opposite. 

Proof. Let f \ 2 be the components of the ground state 
vector-function <p v . According to the Corollary 1, the 
signs of (pi 2 are nx ed. Hence, for the potential V\ 2 of 
fixed sign one may conclude that {ip\ \V\ 2 \ip 2 ) ^ 0. The 
mean value of the Hamiltonian has the form 



+2(rt\V 12 \rt). 



(A17) 



If one replaces <p 2 by —ip 2 and calculates the mean 
value of the Hamiltonian H, the variational principle 
yields immediately (ip\ \ V\ 2 \ip 2 ) < 0. Hence, the signs 
of the product Lp\{x)Lp 2 (x) and V\ 2 (x) are opposite. ■ 

Corollary 3. For the Theorem's conditions the 
ground state is not degenerate. 

[The proof is the same as in one-component case, see 
Ref. Q. Indeed, two vector- functions with components 
which conserve their signs and the sign of the product 
f\(fi 2 cannot be orthogonal. Hence, the ground state is 
not degenerate.] 

Note. If dim(Q,) > 1 one has no rights to insist on 
the validity of the Theorem. However, for the interaction 
potential Vi 2 (x) of a fixed sign one may prove Corollaries 
1 - 3 by similar arguments. 



V. CONCLUSION 

For eigenvalue problems with matrix Hamiltonians the 
structure of nodal points is much more complex than for 
an one-component ID-problem, and the nodal structure 
of multi-component problems is of less importance than 
for one-component ones. Nevertheless, for the interaction 
potential of a fixed sign there are some analogues of the 
oscillation theorems. In particular, for a two-component 
problems for R 1 axis (or a half-axis, or a finite segment) 
one may prove the following statements: 

(1) the ground state is not a degenerate one, and the 
components of the ground-state wavefunction conserve 
their sign, defined by the sign of interaction potential; 

(2) for the state number v (for energies in increasing 
order v — 0,1,2,...) the wavefunction components ip\ 
and ip 2 have the arithmetic mean of the number of nodal 
points (n\ +n 2 )/2 not larger than v. (For a more precise 
form, see Eq. (|A12[> ). 

We suppose it to be sufficient for practical applications. 
Though the condition of a constant sign for Vi 2 seems 
to be too restrictive, one may note, that in practice it 
is sufficient to be sure in a constant sign of V\ 2 (x) for 
the region of localization of the wavefunction components 
only. 
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Relative intensity measurements in the high resolution A 1 T, + ~ 6 3 II — > X 1 E + laser induced 
fluorescence spectra of KCs molecule highlighted a breakdown of the conventional one-dimensional 
oscillation theorem [L.D. Landau and E.M. Lifshitz, Quantum Mechanics, Pergamon, New York, 
1965]. For strongly coupled ^4 1 E + and & 3 II states the number of nodes ua and of the non-adiabatic 
vibrational eigenfunctions ip A and <pl corresponding to the v-th eigenstate differs essentially from 
their adiabatic counterparts. It is found, however, that in general case of two-component states with 
wave functions ip\ and ip\ coupled by the sign-constant potential operator V12 ^ 0: (1) the lowest 
state v = is not degenerate; and (2) the arithmetic mean of the number of nodes ni and 712 of tpi 
and tp? never exceeds the ordering number v of eigenstate: (m + ri2)/2 < v. 
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The conventional oscillation theorem [l], Q states the 
remarkable properties of bound solutions of the one- 
dimensional (ID) Schrddinger equation: (1) none of the 
bound states is a degenerate one and (2) the eigenfunction 
ip v corresponding to the v-th eigenvalue E v vanishes ex- 
actly v times (possesses v internal nodes). Harmonic and 
Morse oscillators, being the basic models for approxima- 
tion of vibrational motion in diatomic molecules, are the 
most famous examples of exact solution of the ID equa- 
tion leading to analytical dependence of eigenvalues E v 
and eigenfunctions ip v on v- values [l], 0] ■ A shape of rovi- 
brational (radial) wavefunction (WF) of a diatomic state 
can be recognized experimentally from the relative inten- 
sity distribution I v (v") in emission spectrum originat- 
ing from a singly excited rovibronic level v'(J'). Nodal 
structure of the upper state WF \<p v | 2 then appears in 
the laser induced fluorescence (LIF) spectra as intensity 
I v (v") oscillation in full progression of vibrational bands 
stretching from (v',v" — 0) to (v',v" — m). This takes 
place in the case of so-called strongly non-diagonal elec- 
tronic transitions, when equilibrium distances r' e and r" 
of the upper and lower states are markedly shifted, see 
Fig. [U Since the vibrational quantum number v' is equal 
to the number of nodes of ip v , the full vibrational LIF 
progression yields absolute vibrational assignment of the 
excited electronic state [3, 0]- The situation becomes 
more subtle when an upper electronic state is subject 
to pronounced perturbations and the adiabatic approxi- 
mation is apparently not valid [5J. In particular, for the 
first excited A 1 Y, + and fo 3 TI states of Rb-containing alkali 
diatomics, strong spin-orbit (SO) interaction matrix ele- 
ment £j£b l ea ds to apparent disorder in vibrational spac- 
ing observed in 0] for Rb2 and for NaRb. Though in 
such a situation vibrational numbering becomes a compli- 
cated task, unambiguous w'-assignment of NaRb A-state 
was achieved in 0,111 by analyzing LIF intensity distri- 
bution I v (v") in long A ~ b — > X progressions. This 
became possible since the number of maxima in observed 
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FIG. 1: (color online) Schema of potential energy curves for 
isolated ground X 1 E + and deperturbed excited ^4 J E + and 
& 3 II states of KCs molecule. Inset represents the relevant 



spin-orbit coupling function £%b(r). 



intensity I v 



precisely follows the conventional ID os- 



cillation theorem. In case of alkali diatomics containing 
heavy Cs atom the situation becomes more complicated 
0-11 1 because of the adiabaticity parameter £™ b /y/LJAWb 
is close to one for the A ~ b complex, uja, being har- 
monic frequencies. A significant mixing of adiabatic (or 
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FIG. 2: (color online) Relative intensity distributions in the 
A ~ b(E v '-f) -»■ X(v'x,J" = J' ± 1) LIF spectra of KCs 
originated from the J' = 50, 109 and 128 levels with predom- 
inant fraction of the = state. Numbers over bars denote 
vibrational quantum number of the ground X-state v' x ■ The 
inset zooms additional maximum for J' = 50. 



FIG. 3: (color online) The two-component non-adiabatic 
wavefunctions tp v A (r) and ip% (r) calculated for the rovibronic 
levels of the KCs A £ + ~ & 3 FI complex presented in Fig[2] 
Xa (r) is diabatic WF of the deperturbed yl-state for — 0. 
Insets enlarge a nodal structure of the non-adiabatic (p\ (r) 
wavefunctions in the region r £ [3.5,4.5] A. 



diabatic) vibrational WFs [5J should inevitably affect the 
shape of the non-adiabatic WFs; indeed, peculiarities in 
the A ~ b — > X intensity distribution have been already 
observed for NaRb [l2| and NaCs molecules. 

In present study we report on the breakdown of os- 
cillation theorem in the A ~ b — » X LIF of the KCs 
molecule unambiguously observed even for the lowest vi- 
brational level of the A-state, and suggest, basing on the 
model study, a non-adiabatic analogue of oscillation the- 
orem which is valid for two-component states coupled 
by the sign-constant potential operator V12 7^ 0. It 
should be noted that polar alkali diatomics are of fun- 
damental interest due to recent progress in the produc- 
tion and trapping of ultracold species 13j, ll4[ leading to 
their possible applications in quantum information de- 
vices controlled by an external electric field [15j . Favor- 
ably for the above applications, the singlet-triplet A ~ b 
complex can provide with intermediate levels for effi- 
cient absorption-emission cycles into "absolute" ground 
state a 3 E+(Ef,J" = 0) ->■ A ~ b(E v ';J' = 1) 

= 0, J" = 0) for formation of stable ultracold 



molecules [161 ]. 

The experimental setup and details of measurements 
can be found elsewhere [13, E3] • Briefly, KCs molecules 
were produced in a stainless steel heat pipe and kept 
at about 280°C temperature. The scanned diode lasers 
with 850 nm, 980 nm and 1020 nm laser diodes were 
employed for A ~ b <— X excitation. Backward LIF 
was collected on the input aperture of the Fourier trans- 
form spectrometer Bruker IFS 125HR with the resolution 
0.03 cm -1 leading to the uncertainty of the line posi- 
tions about 0.003 cm -1 . Relative intensity distributions 
were determined with about 5% uncertainty taking into 
account the spectral sensitivity of the InGaAs detector. 
Fig. [5] presents A ~ b — > X LIF intensities from the 
levels with J' — 50, 109 and 128 with respective term 
values 10188.33, 10401.79 and 10508.09 cm" 1 related to 
the minimum of the ground state potential energy curve 
[l7j |. These levels belong to the lowest vibrational level 
= of the diabatic A-state with predominant sin- 
glet fraction P\ = {tp\\<p\) r being 88.6%, 81.0% and 
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79.6%, respectively. This was determined [18[ by a de- 
perturbation analysis of the experimental term values of 
the KCs A ~ b complex in the framework of the inverted 
channel-coupling approach by means of the 4x4 model 
Hamiltonian constructed on Hund's coupling case (a) ba- 
sis functions similar to [Io|, 12 1 , allowing constructing di- 
abatic A 1 Y> + and fe 3 IIn potentials and relevant SO cou- 
pling function (see Fig. [IJ and to reproduce term values 
with experimental accuracy of about 0.005 cm" 1 . In- 
deed, a single maximum of intensity distribution in the 
LIF spectrum originating from the level with J' = 128 
confirms the energy based v A = assignment. However, 
manifestation of a smaller second maximum observed in 
the region of 8900 cm -1 for the levels J' = 50 and 109 
unambiguously highlights the presence of at least one ad- 
ditional oscillation in a relevant upper state WF, thus, at 
least for the J' = 50 and 109 levels there is an apparent 
contradiction between the observed intensity distribution 
and the vibrational assignment based on the deperturba- 
tion arguments. 

To elucidate such a discrepancy the A ~ b(E" ; J') — > 
X(v";J" = J' ± 1) intensities I\^b-^x( v ") were simu- 
lated as 



TV 

J A~b^X 



W) 



VA^b^rX 



u A~b^,X 
jpv' 



W A \dAx\x v x. 



(1) 



&x > 



where <Iax ( r ) is the spin- allowed — transi- 

tion dipole moment, VA~b-yX is the transition wave num- 
ber, E A ^ b is the energy of the rovibronic level v'(J') of 
the A ~ b complex, ip A (r) is the singlet A-state frac- 
tion of non-adiabatic WF, while E x and Xx ( r ) are the 
adiabatic energy and WF of the X-state, respectively. 
As follows from Eq. (1), the intensities provide infor- 
mation essentially about the singlet A-state WF since 
the b 3 H — > X 1 Yi + transition is strictly forbidden, i.e. 
dbx = 0. To obtain the required non-adiabatic WF the 
simplified two-channel form of Hamiltonian [19j 



hi Va 
V 2 i h 2 



= E v 



(2) 



was used. Here, v € [0, N] is the ordering number of the 
eigenvalue E v , hi are the ordinary ID radial Hamiltoni- 
ans 



1 ~ 2fi dr 2 



Ui(r) 



ie[l,2] (3) 



coupled by Hermitian operator Vi2(r), /i is the reduced 
mass, Ui are diabatic rotationless potentials while bound 
state eigenfunctions ip\ obey the conventional bound- 
ary ifii(0) = Q\<p1(r — > oo) —r and normalization 
PI + P% = 1 conditions, where P? = (<^|^)r- In- 
serting into Eq. ([2]) the deperturbed diabatic poten- 
tials Ui = UaW), U2 = Ub{r) and SO coupling function 
V12 = V2£,Ab ( r ) UM > see Fig- E the non-adiabatic eigen- 
values E v = E\^ h and eigenfunctions <~p\ = ip A (r) and 
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FIG. 4: (color online) Number of nodes tia, of two- 
component vibrational wavefunctions of the KCs A ~ b com- 
plex ip A (r) and ip% (r) counted as a function of ordering num- 
ber v G [0, 1, ••• ,N] and perturbation parameter A. The 
straight line corresponds to the limit case tia — nt — v' . 



ip\ = ipl (r) of the A ~ b complex were obtained. The 
obtained WFs are presented in Fig.[3j along with the dia- 
batic A-state WFs Xa(t) corresponding to v A = level. 
The present deperturbation model accounts only for the 
SO interaction between ^4 1 S + state and 6 3 Hq = o compo- 
nent of the triplet state as applied in @, Q • Inserting in 
Eq. ([J) the energies E A „ b and WFs tp v A along with ab 
initio oIax [H| and adiabatic eigenvalues E x and eigen- 
functions Xx °f the ground state from [T3] yielded the in- 
tensity distributions shown in Fig. [3] Remarkable agree- 
ment of the predicted intensities with their experimental 
counterparts justifies the applicability of the simplified 
2x2 deperturbation model. As can be seen from Fig. [3[ 
for J' = 128 the shape of non-adiabatic WF ip A (r) in the 
classically allowed region of motion r € [4.5, 5.5] A is sim- 
ilar to a diabatic one, thus confirming the v A = assign- 
ment, see Fig. [21 additional WF minimum at r = 5.33 A 
has a very small amplitude of 0.006 and can not be seen in 
the scale of Fig. [3] In contrast, WFs ip v A {r) for J' = 109 
and 50 have a pronounced additional minimum and a re- 
spective node located at r ~ 5.1 A which is responsible 
for the second maximum in the respective intensity distri- 
butions in Fig. [5] thus explaining the observed breaking 
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of the oscillation theorem. Careful analysis (see insets 
in Fig. [3]) reveals that the full number of nodes ua of 
the non-adiabatic WF ip\ for J' = 50, 109 and 128 lev- 
els possessing respective ordering numbers of eigenvalues 
v' = 20, 19 and 19, is equal to ua = 14, 13 and 13 respec- 
tively. However, most "additional" nodes of ip\ appar- 
ently could not be observed in the LIF spectra since they 
are located in the classically forbidden region of motion 
of the deperturbed A-state where the amplitude of WF is 
very small, see inserts in Fig. [3l The non-adiabatic tpl in 
Fig. |3] for J' = 50, 109 and 128 possesses the respective 
number of nodes rib — 16, 15 and 15 which is close to the 
diabatic 6-state vibrational number v* % — 16, 15 and 14 
as discovered from rigorous deperturbation analysis [3] ■ 
To find an analogue of oscillation theorem for a two- 
component system y\ and ^ it is necessary to establish 
a connection between the respective number of nodes ni , 
ri2 and the ordering number of the eigenstate v. For this 
purpose it is useful to trace such a connection for the 
KCs A ~ b system under study at different strength of 
the SO interaction V12 = A x V2^ A °b which is varied by 
a factor A. To get rid of the influence of a particular 
J' value determining the effective potential, the simu- 
lation of ip\ and (fil yielding the respective number of 
nodes tla and rib was performed for the rotationalless 
case J' = 0. The resulting u'-dependencies of ua and 
rib values are presented on Fig. |4] for different A values. 
It can be seen that, at fixed A, for all v' below some 
particular v max , the number of nodes coincides with v': 
n a = rib = v' < v max , where v max monotonically in- 
creases with A being v max — 19, 33 and 57 for A = 1, 
5 and 10, respectively. The smallest number v max = 17 
corresponds to vanishing perturbation A — > 0. The cases 
riA = rib < v 1 and ua 7^ rib *C v' take place for high 
levels v' > v max - The quasi-periodic alterations of ua 
and rib observed at v' > v max is attributed to acciden- 
tal degeneration of deperturbed levels of diabatic A and 
b states [3]. The most important feature is that in all 
cases the arithmetic mean of the number of nodes ua and 
rib never exceeds the ordering number v'. This general 
property allowed us to formulate a non-adiabatic ana- 
logue of the oscillation theorem valid for two-component 
states mutually perturbed by a sign-constant potential 
operator V\?,{r) ^ 0: (1) the lowest bound v = state is 
not a degenerate one and the respective two-component 
WFs ip\ =0 , <^2 = ° do n °t have nodes; (2) the number of 
nodes n\ and ri2 of the wavefunctions ip\ and ip^ corre- 
sponding to the v-th non-adiabatic eigenvalue obeys the 
inequality n\ + ri2 < 2u. Rigorous mathematical proof 
of the theorem is obtained basing on minimax principle 
[20I ] and is presented in the Appendix. In practical imple- 
mentation of the theorem a constant sign of the coupling 
operator is important only for the limited region where 
the WFs are localized. In particular, it is easy to verify 
that the theorem is valid for the rovibronic J' = 50, 109 
and 128 levels of the KCs A ~ b complex studied above 



since ua + rib < 2u' (see Fig. [3]). 

Thus, the spin-orbit coupling effect can affect the nodal 
structure of a multi-component diatomic vibrational WF 
leading to a breakdown of the conventional ID oscillation 
theorem. Nevertheless, for two-component states coupled 
by the sign-constant potential operator non-adiabatic 
analogue of the oscillation theorem has been formulated 
and proved. 
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APPENDIX: Proof of Oscillation Theorem for 
Nonadiabatic Diatomic States 

I. NOTATIONS AND GENERAL NOTES 

Here we consider the two-component Schrodinger 
equation (see Eq.(2) of the Main Text) for the 2x2 ma- 
trix Hamiltonian H: 



hi V12 
V21 h 2 



E v 



£V- (Al) 



The components h\p of H in Eq. (|Al[) have the tra- 
ditional form hk = — ^7^2 + Uk(x) with some effec- 
tive mass fi and real- valued potentials Uk (here the ro- 
tational part of energy is included in Uk)- The inter- 
action operators are supposed to be some real-valued 
potentials V2i(x) — Vi2(x), and the energy levels E v 
are enumerated in increasing order, taking degeneracy 
into account, by traditional vibrational quantum number 
v = 0, 1, 2, • • • . The matrix form of Eq. (|Alj) is equivalent 
to the following system of equations: 



21^1 + ^2f2 = & W 



(A2) 



For molecular problems one may suppose the compo- 
nents (pi 2 ( x ) °f a vector- function <p v to be real- valued 
functions defined for some region f2, being the same for 
both components. It is supposed here O = [a, 0\ C M 1 
for some finite (or infinite) a and (3. When one of the 
ends of the interval [a, /?], for example, a, is finite, one 
supposes that both components of solution ip v obey the 
Dirichlet boundary condition at this point: y>\ 2 (pi) = 0. 



5 



In any case only bound states are considered and solution 
ip v is supposed to be normalized by condition 



(<pm = (<Pi\<p?)+(<Pt\<pZ) = i. 



(A3) 



It is worth to note that for invertible potential V12 (x) one 
may use the system Eq. (IA2p to find (pj as a solution of 
the 4-th order differential equation (similar relation holds 
for ipl): 



[(h 2 - E^V^ihx - E v ) + V l2 ]<pX = 0. 



(A4) 



In particular, as a result of a well known property of 
the Cauchy problem, the following Lemma is useful for 
analysis of Eq. (|A2p : 

Lemma. For invertible potential Vi2(x) the compo- 
nents of any nontrivial solution ip v of the problem (IA1|) 
cannot be identically zero into some open subregion of 

n. 

For the ID one-component problems it is impossible 
that the wave function and its gradient vanish simul- 
taneously at some point. In particular, near the nodal 
points the wavefunction changes the sign. For the multi- 
component problems situation is not so simple. It is pos- 
sible that at some point some component of the wave 
function only touches the axes. When the interaction 
potential does not vanish at this point, it follows from 
Eq. (|A2[) that both components (p\ 2 (x) vanish along with 
their second derivatives. This special situation is not sta- 
ble with respect to small perturbations and we ignore it 
here. That is, we consider as the nodal points only the 
ones for which functions change the sign. Similar issues 
for one component problems for 51 = M. n (n > 1) enforce 
to analyze the regions of a constant sign for wave func- 
tions [2J. 

II. SOME STATEMENTS 

For one-component one-dimensional problems it is well 
known that the vibrational state number v coincides with 
the number of nodal points (see Ch.XIII.3 in [20| or 
Ch.VI, §6 in [Hj]). For matrix equation (|Al|) the problem 
is much more difficult. 

Let <p v be an eigenvector of the two-component ma- 
trix Hamiltonian H with energy E v and components <~p\ 
and (f2- Let us analyze the regions of a constant sign for 
functions <^i 2) that is, the intervals where the function 
has a definite sign and changes it at the ends. The num- 
ber Kj of such regions corresponds to rij = Kj—1 nodal 
points of (p'j (j = 1,2). Let us suppose that ip\ has K\ 

regions fl^. (j = 1,2, ...,K±) of a constant sign. One 

may define for each of regions Q,j the function ip^ that 

equals to tp\ on Clj and equals to zero out of flj . Let 

(2) 

us suppose that (p\ has K 2 regions Q k ' [k = 1,2,..., K 2 ) 

(2) 

of a constant sign and ip k are the functions similar to 



the ones for the first component of (p v . Hence ip v defines 
the linear space JC of (K\ + K2) -dimensional vectors C 
defined by any set of K\ coefficients Cj and set of if 2 co- 
efficients dk and corresponding vector-functions \ with 
components defined by the relation 

K x K 2 



Xi 



E 



Cj <p 



(1) 



X'2 



= E rf 



(2). 



(A5) 



fc=i 



Note that xi,2 and their gradients have to be the square 
integrable functions on f2. 

Here we will study some properties of the matrix M, 
defined for the above mentioned space K by the vector- 
function x in the following way: 

(C,MC) = - 4) 2 ^f|F 12 #i 2) ). (A6) 

J, As 

The described notations are used in all statements of this 
section. We use also the symbol [n/2] for the integer part 
of the n/2 value. The following relation is important for 
our discussions: 

" +1" 



(A7) 



Our further discussion is based on the simple expression 
of energy functional E(x) associated with the Hamilto- 
nian H, for the vector-function x with components of the 
form (|A5]l . 

Statement 1. If ip v is a solution of Eq. (|Al[) . then for 
any vector C £ K. one may write 

E(x) = EV(x\x) - £( Cj - duf^plV^)- (A8) 

j,k 

Proof. The use of a standard form for the energy func- 
tional (with kinetic energy expressed as squared gradi- 
ent of wavefunction) and integration by parts enables to 
write, due to Eq.(|A5p. 



K-2 



^)=E c ?^ 1) i /i ^f ) )+E4(4 2) i^i^ 2) ) 



fc=i 



+2^c j d k ( l pf ) \V 12 W k 



(2)\ 



(A9) 



where integration in each of integral in the sums is done 



over the regions ^l)f' or ty 1 ' f] f2^, respectively. 



)(2 



(1), 



)(2) 



The use of Eq. (IA2[) gives immediately 



E(x) = E v 



K-2 



2/„(2) v ( 2)) 



fe=l 



E 



K 2 



fc=i 



(2)\ 
k I 

( 1 )|„(2)\ 



(A10) 
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This is exactly expression 

Statement 2. There are no more than (K\ + K 2 — 1) 
regions Qj |~| fi^ C f2, where the product (p\^p 2 differs 
from zero and conserves the sign. 

Proof. Let us suppose that there exist Kj regions of a 
constant sign for <p° . There are no more than nj = Kj — 1 
internal points of f2 where <pj changes the sign. Hence, 
there are no more than ni+n2 points where f\(f2 changes 
the sign. That is, there exist no more than (Ki — 1) + 
(K.2 — 1) + 1 = (if i + K-2 — 1) non-empty regions of the 
type Clj P| Ctfjp where ¥>3V2 conserves the sign. ■ 

Statement 3. Let us suppose that the sign of the 
potential Vn{x) is fixed. Then the matrix M defined 
by Eq. (|A6l) has no less than [ " R " 1+ 2 ^ 2+1 ] non-negative 
eigenvalues. 

Proof. According to definition, the quadratic form 
(C,MC) is determined by a system of matrices of 
rank 1, each of which has only one nonzero eigenvalue 



and, hence, 



2(tpy\Vi2\tp^} ■ According to Statement 2, among those 
eigenvalues there are no more than (K\ + K2 — 1) non- 
zero values; half of them is negative, while another half is 
positive. In any case the largest number of strictly neg- 
ative integrals is no more than (see Eq. (|A7p \ Kl % I<2 ] ■ 
According to the minimax principle (see, e.g. Ch.VI. of 
[22| or Sect. XIII. 1 of [111), this means that no more than 
[ 2 +1 ] eigenvalues of matrix M are strictly negative. 
Hence, the remaining eigenvalues of M are non- negative 
and their number is no less than 



(K 1 + K 2 ) - 



~K X + K 2 ~ 




~K X +K2 + I 


2 




2 



(All) 



III. THE NODAL POINTS IN 
TWO-COMPONENT PROBLEM 

For ID problems the wavefunction of the state with 
vibrational number v = 0,1,2,... has exactly v nodes 
H[2lj]. When dim(fl) > 1 for one-component problem 
one may only prove that, for a given state, the number of 
regions of a definite sign is not too large [2J . For the two- 
component system considered here one may prove some 
analogue of these statements. The proof is similar to 
the one in the one-component case, but it requires some 
special modifications. 

Theorem. For the stationary state E v , tp v of the 
problem (|A1|) with the interaction potential V12 of a def- 
inite sign the number of points n\ and n 2 , where the 
respective components Lp\ and tpi\ change their signs, the 
following inequalities hold: 



ni + n 2 + 1 



(A12) 



m + n 2 



< v. 



(A13) 



Proof. Let us consider the space K, associated with re- 
gions of constant sign for components of the wavefunction 
ip v , as was described in Section [TIJ Each vector C E JC 
may define some vector- function \ (see Eq. (|A5| ). due to 
Statement 1 the corresponding energy functional value 
can be written as 



£(X) = ^<X|X>-(C,MC). 

According to Statement 3 (note that, Kj 
matrix M has no less than 



T 



(A14) 
1), the 



K x + K 2 - 


rl" 




ni + n 2 + 3 


2 




2 



n\ + n 2 + 1 



(A15) 



non-negative eigenvalues, each of them being associated 
with some vector C and corresponding vector-function \- 
Let us suppose T > v + 2. Then one may find at least 
two independent normalized vectors x^ 1 ' 2 -* (or, equiva- 
lently, C^ 1,2 ' € /C being orthogonal to v exact solutions 



of Eq.jATj) with lowest energies. For 
these vectors (C^MC^)) > (j = 1,2) and, hence, 
E(x( 1,2 )) < E v . It follows from the variational princi- 
ple (see, e.g. Ch.VI. of 22] or Sect.XIII.l of [20]) that 
both functions an d are solutions of Eq. (|Al[) with 
energy E v and one may construct, as a linear combina- 
tion of these vector-functions, a non-zero solution x 01 
Schrodinger equation, for which, at least for one of re- 
gions Clj or fijj. , the corresponding component of the 
solution (xi or X2) identically equals to zero. This is 
however impossible, according to the Lemma of Section 
[TJ Hence, T < v + 2, or, equivalently, T < v + 1, and one 
may write 



ni + n 2 



< 



"1 



n 2 



< v. 



(A16) 



The use of Eq. (IA7[) proves the statement. ■ 

Corollary 1. For the Theorem's conditions, the com- 
ponents of the function (p" =0 conserve their signs into il. 
Indeed, for v — inequality (IA13[) means n\ = 112 = 0. 

Corollary 2. For the Theorem's conditions for the 
ground state and V\ 2 {x) < the signs of (p\ and tp^ 
coincide, while for Vi 2 (x) > the signs are opposite. 

Proof. Let ip\ 2 be the components of the ground state 
vector- function <p v . According to the Corollary 1, the 
signs of (pi 2 are fixed. Hence, for the potential Vi 2 of 
fixed sign one may conclude that (ip\ \V12\1p2) 7^ 0. The 
mean value of the Hamiltonian has the form 



+ (<f%\ha\<f%) 



(A17) 
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If one replaces ip^ by —ip^ an d calculates the mean value 
of the Hamiltonian H, the variational principle yields 
immediately {tp\ \V12\W2) < 0- Hence, the signs of (fi\,<fi2 
and V12 are opposite. I 

Corollary 3. For the Theorem's conditions the 
ground state is not degenerate. Indeed, two vector- 
functions with components which conserve their signs 
and the sign of the product f\^>2 cannot be orthogonal. 
Hence, the ground state is not degenerate. See also Ref. 
0. It is clear, that (n\ + ri2)/2 is the arithmetic mean of 
number of nodes for components of the two-component 
solution. Hence, the Theorem is an analogue to the usual 
one-component statement. 

Note. If dim(n) > 1 one has no rights to insist on 
the validity of the Theorem. However, for the interaction 
potential Vi2{x) of a fixed sign one may prove Corollaries 
1 - 3 by similar arguments. 

IV. CONCLUSION 

For eigenvalue problems with matrix Hamiltonians the 
structure of nodal points is much more complex than for 
a one-component ID problem, and the nodal structure 
of multi-component problems is of less importance than 
for one-component ones. Nevertheless, for the interaction 
potential of a fixed sign there are some analogues of the 
oscillation theorems. In particular, for a two-component 
problems for M 1 axes (or a half-axes) one may prove the 
following statements: (1) the ground state is not a de- 
generate one, and the components of the ground-state 
wavefunction conserve their sign, defined by the sign of 
interaction potential; (2) for the state number v (for en- 
ergies in increasing order v — 0,1,2,...) the wavefunction 
components <p\ and have the arithmetic mean of the 
number of nodal points (n\ + n2)/2 not larger than v. 
We suppose it to be sufficient for practical applications. 
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